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We present systematic temperature-quench Monte Carlo simulations on discrete-strain pseudospin 
model Hamiltonians to study microstructural evolutions in 2D ferroelastic transitions with two- 
component vector order parameters {Nqp = 2). The zero value pseudospin is the single high- 
temperature phase while the low-temperature phase has Ny variants. Thus the number of nonzero 
values of pseudospin are triangle-to-centered rectangle (W = 3), square-to-oblique {Ny = 4) and 
triangle-to-oblique {Ny — 6). The model Hamiltonians contain a transition-specific Landau energy 
term, a domain wall cost or Ginzburg term, and power-law anisotropic interaction potential, induced 
from a strain compatibility condition. On quenching below a transition temperature, we find be¬ 
haviour similar to the previously studied square-to-rectangle transition {Nqp — 1, Ny — 2), showing 
that the rich behaviour found, is generic. Thus we find for two-component order parameters, that 
the same Hamiltonian can describe both athermal and isothermal martensite regimes for different 
material parameters. The athermal/isothermal/austenite parameter regimes and temperature-time- 
transformation diagrams are understood, as previously, through parametrization of effective-droplet 
energies. In the athermal regime, we find rapid conversions below a spinodal like temperature and 
austenite-martensite conversion delays above it, as in the experiment. The delays show early incu¬ 
bation behaviour, and at the transition to austenite the delay times have Vogel-Fulcher divergences 
and are insensitive to Hamiltonian energy scales, suggesting that entropy barriers are dominant. 

PACS numbers: 64.60.De, 81.30.Kf, 64.70.K-, 05.70.Ln 


I. INTRODUCTION 

Steels and shape memory alloys are martensitic ma¬ 
terials that undergo diffusionless, first-order phase trans¬ 
formation from high-temperature parent ’austenite’ unit¬ 
cell to low-temperature product ’martensite’ unit-cells 
(or variants) on cooling or under external stres^. A sub¬ 
set of physical strains are the order parameter (O P)- As 
martensitic materials have many applications^, much 
work has been done to understand domain-wall mi- 
crostructures and their underlying kinetics. According 
to traditional classificatioiP, martensites are classified 
as athermal, with rapid milli-second austenite-martensite 
conversions on cooling below a martensite start tem¬ 
perature and no conversions above it; and isothermal, 
which can have slow or delayed conversions in minutes or 
hours. But, experiments on athermal martensitic materi¬ 
als have found delayed conversions above the martensite 
start temperature, where only austenite should existP. 
Computer simulations of martensitic models could give 
insights into the classification of martensites and the un¬ 
expected delayed-conversions in athermal martensites. 

Continuous variable nonlinear free energies are mini¬ 
mized in displacement, phase field, and strain using relax- 
ational dynamics, Monte Carlo (MC) and Molecular dy¬ 
namics simulation^^H^ and the obtained microstructures 
are consistent with experimenfP^, but that can need ex¬ 
tensive computer time. More economic discrete-strain 
clock-like model Hamiltonian^^ are systematically de¬ 
rived from continuous strain free energies for different 
ferroelastic transitions in 2- and 3-spatial dimensions 
(2D & 3D). Power-law anisotropic interaction potentials. 


which ar ise fr om the no-defect St.Venant compatibility 
conditiorPm^, and orient strain domain walls, have their 
counterparts induced in pseudo spin Hamiltonians. The 
microstructures generated from these strain-pseudospin 
models using local mean-field approximatioiP^ are in 
good agreement with the continuous variable modelJ^^^ 
and experiments^. 

Systematic temperature-quench MC simulations were 
performed on the simplest scalar-OP, 3-state pseudospin 
Hamiltonian for square-to-rectangle (SR) transit iorP^l 
and showed both rapid conversions below a spinodal- 
like temperature and incubation-delays above it, as in 
experiments^ on athermal martensitic materials. The 
conversion-time delays found to have Vogel-Fulcher di¬ 
vergences, which are insensitive to Hamiltonian energy 
scales and log-normal distributions, suggesting the dom¬ 
inant role of entropy barriers. An athermal/isothermal 
martensites regime diagram is predicted in material- 
parameters; crossover temperatures and domain-wall 
phases in Temperature-Time-Transformation (TTT) di¬ 
agrams are understood through parametrization of tex¬ 
tures by surrogate droplet energies; and role of power- 
law potentials are shown to be important for textures 
and incubation^. The central question is: Are such 
conversion-delays in the athermal martensite regime, spe¬ 
cific to the scalar-OP transition, or are they generic, ap¬ 
pearing in vector-OP transitions ? 

In this paper, we show that the athermal martensite 
regime conversion-delays in the scalar-OP {Nqp = 1) 
SR transition are generic in three vector-OP {Nqp = 
2) ferroelastic transitions: triangle-centered rectangle 
{Ny = 3); square-oblique {Ny = 4); triangle-oblique 
{Ny = 6). Under systematic MC temperature quenches. 
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we find isothermal parameter regime with slow or delayed 
conversions and athermal parameter regime that has 
rapid conversions below a temperature and incubation- 
delays above it, as in experimeniP and scalar-OP SR 
transitiorP^. The athermal regime conversion-time de¬ 
lays have Vogel-Fulcher divergences, which are insensi¬ 
tive to Hamiltonian energy scales and log-normal dis¬ 
tributions. The athermal/isothermal/austenite regime 
diagrams are obtained in material parameters. The 
crossover temperatures and domain-wall phases in the 
TTT diagram are understood through the parametriza- 
tion of textures. Microstructures obtained in these tran¬ 
sitions are in good agreement with continuous-variable 
simulation^^H^ and experimenfP^. We finally show the 
importance of power-law interaction potentials in the in¬ 
cubation behaviour, and microstructures. 

The paper is organised as follows. In Section 2 , we 
outline derivations of the vector-OP strain-pseudospin 
Hamiltonians. In Section 3, we present the athermal/ 
isothermal martensite regimes and crossover in mate¬ 
rial parameters. In Section 4, we focus on the athermal 
martensite regime and present conversion-delay kinetics, 
parametrization of domain-wall phases in TTT diagram 
by effective droplet energies, and conversion incubation 
textures. In Section 5, we present kinetics in the absence 
of the power-law anisotropic interactions that shows de¬ 
lays without incubation, and Section 6 is a summary. 


II. STRAIN-PSEUDOSPIN HAMILTONIANS 

In this Section, we state for completeness, the vector- 
OP strain-pseudospin model Hamiltonian^, that were 
systematically derived from scaled continuous-strain 
free-energies^ for triangle-to-centered rectangle (TCR), 
square-to-oblique (SO) and triangle-to-oblique (TO) fer- 
roelastic structural transitions. 

In 2D, structural transitions have d{d + l)/2 = 3 or 
three distinct physical strains, the compressional (ei), 
deviatoric ( 02 ) and shear ( 03 ) strains. Of these, ( 02 , 63 ) 
are OP {Nqp = 2 ) and the ei is non-OP {n = 1 ) strains. 
The scaled free energy has a Landau term Fp ; a Ginzburg 
term, quadratic in the OP gradients Fq; and a seemingly 
innocuous term, quadratic in the non-OP strains 
that turns out to generate crucial power-law anisotropic 
interactions between the OP strain^^^. Thus 

F = Eo[Fl + Fq + Fnon] (2-1) 

Here Fq is an elastic energy per unit cell. The transition 
specific Landau term Fp has {Ny + 1 ) degenerate energy 
minima at the first-order transition as shown in Figure 
1. The high-temperature austenite minima is allowed at 
all temperatures as its existence has to be determined 
dynamically, and Ny minima are the low-temperature 
martensite variants. The pseudospin derivation results 
of Ref 11 are restated here, for completeness. 






FIG. 1. Landau free energy minima and strain-pseudo spin 
eloek veetors: (a) Contours of Landau free energy Fl showing 
single austenite minima at center and Ny martensite minima 
at corners, in a plot of es versus 62 and (b) corresponding 
clock model minima at pseudospin values for TCR (left), SO 
(middle) and TO (right) transitions. 

The scaled Landau free energy Fp for TCR transitiorl^ 
Fl{^ = ^(r - 1)^ + [^ - 2{el - Segel) + (e^)^], (2.2) 

f 

has an austenite minima at ( 02 , 03 ) = (0, 0), and Ny = 3 
martensite minima at which ( 02 , 03 ) = (cos 0, sin 0) for 
0 = 0,^,^. Here^ =4 + el and r = (T-Tc)/(To- 
Tc) is the scaled temperature; Tq is the first-order Lan¬ 
dau transition temperature and Ty is metastable austen¬ 
ite spinodal temperature. 

The scaled Landau free energy Fp for SO transitiorP^ 

= + 4e® - c'ielel, (2.3) 

r 

also has an austenite minima at (02,03) = (0,0), 

and Ny = 4 martensite minima at which (02,03) = 

(cos 0, sin0) for 0 = J, ^ and ^ with material de¬ 
pendent elastic constant C^. 

The scaled Landau free energy for TO transitioiPIl is 

Fpe) = (2.4) 

r 

where Cq is a material dependent parameter. The Lan¬ 
dau polynomial has an austenite minima at ( 02 , 03 ) = 

(0, 0), and Ny = 6 martensite minima at which (02, 03) = 

(cos(/),sine/)) for (/> = 0, f, ^ and 

The domain-wall cost Ginzburg term Fq is, 

^G(Ve-)=^2^(Vey (2.5) 

r 

The non-OP term is harmonicPH, with stiffness Ai, 

F„.„(ea) = ^Fe2 = ^T|e^|2 ( 2 , 6 ) 

r k 
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and is minimized subject to St.Venant compatibility con¬ 
straint for physical strain^, 

A^ei - {Al - Al)e 2 - 2A^Aye^ = 0; (2.7a) 

with gradient terms as difference operators V ^ A for 
sites r on a computational grid. In Fourier space 
K^(k) = 2sin(/c^/2) and so 

Oiei + 02^2 + Os^s = 0 (2.76) 

where the coefficients for square lattice are Oi = 

O 2 = — Ky)^ and O 3 = 2KxKy; for triangle lat¬ 
tice, Oi = O 2 = {Kl-Kl), and O 3 = ‘^K^Ky. 

Here, = {K‘^ + ^y)' Minimization of non-OP strain 
generates power-law anisotropic interactions between OP 
strains, by inserting a direct solution ei = —( 02^2 + 
0363)101 for ^ 7 ^ 0 into ( 2 . 6 ), 

Fcompat{.e 2 ,ez) = ^ y] et{k) U^'ik) e*,{k) (2.8) 

£/=2,3-k 

where U 22 {k) = y{P 2 /Oxf , I/23(^) = y{ 0203 l 0 rf , 
U 33 {k) = y{ 03 / 0 i)‘^ and u = (1 — (5^^). Figure 2 shows 
the power-law potentials as relief plots in Fourier 
space and contours in coordinate space. 



FIG. 2. Power-law anisotropic potentials: Relief plots of ker¬ 
nels U 22 {k), U 23 {k), and Ussik) in Fourier space {top row) and 
corresponding contours in coordinate space {bottom row). 


The continuous-strain OP e = { 62 , 63 ) is discretizecff^ 
by choosing its values only at the Ny + 1 Landau minima, 

e(r) = |e|(cos0, sin0) s{r)S{r). (2.9) 

The Landau ternffH becomes. 


where = r - 1 + {s - 1)^,5^(t) = |{1 -h 

for TOR, and gp = r — 1 -h (^^ — 1)^,5 ^(t) = |{1 + 

yT^^'SrTd} for SO and TO transitions. 

The square gradient Ginzburg term becomes, 

Fg(V 5) =C'^e'(V5)2 = ^2^772g-2|5(fc)|2 (2.11) 

r k 

The discrete-strain^ _ pseudospin clock-zero model 

Hamiltonian is derived^ by substituting (2.9) into the 
total free energy (2.1), 

PH{S) = PF{e^sS) (2.12) 

The Hamiltonian in coordinate space is 

PH = ^[Y^{gPr)S^{r)+e{^F} 

f 

+ S T! ^Uu'{r-r)Si{r)Se{r)\, (2.13a) 
r,r' £/=2,3 

where Dq = 2 EQe‘^. It is diagonal in Fourier space, 

= E Qo,ee'ik)Se{k)S}{h (2.136) 

k £/=2,3 

= Do[{ 9 LiT)+eK^}Se/ +Fuu,ik)], (2.13c) 

and is a clock-zero model Hamiltonian with single austen¬ 
ite S = {S 2 , S 3 ) = (0, 0) and Ny martensite variants: 

s= (l,0),(--,±^);(±-,±-);(±l,0),(±-,±^) (2.14) 

for TCR (A„ + l = 4), SO {N^+1 = 5), and TO (JV^ + l = 

7) transitions respectively. 

MC temperature-quench simulations are carried out 
systematical!}!^ on a square lattice in 2D. At t = 0, we 
consider 2% of sites randomly with Ny strain-pseudospin 
martensite values in austenite. The seeds are quenched 
below the Landau transition r << 1 and held for t < tp 
MC sweeps (MCS). Metropolis algorithm is used for ac¬ 
ceptance of energy changes, that are calculated through 
Fast Fourier transforms. In each MC sweep, we visit all 
N = L X L sites randomly, but only once. Simulation 
parameters are L = 64, Tq = 1; Tc/Tq = 0.6, 0.7,0.8,0.9, 
e = l;Ai = 1,4,10;2 Ai/A3 = 1; F;o = 3,4,5,6;4 < 
10,000 sweeps, and conversion times are averaged over 
Nyuns = 100 runs. 


III. ATHERMAL AND ISOTHERMAL 
PARAMETER REGIMES 


Hl{S) = e^J29L{T)S^{f) =e^J29L{T)\S{k)\^ (2.10) 

r k 


On quenching 2% of martensite seeds to a temperature 
r(T) < t(To), we defin^i^ martensite conversion fraction 
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FIG. 3. Crossover behaviour of athermal/isothermal marten¬ 
site eonversions: Left: Conversion times tm vs T/Tq show¬ 
ing crossover from athermal (fast) to intermediate (slow) for 
fixed = 3,Ai = 4 and Tc = 0.9, 0.8, 0.7. Right: Conver¬ 
sions showing crossover from intermediate (slow) to athermal 
(fast) for fixed Eq = 3,Tc = 0.7 and Ai = 4,2,1. Top row: 
TCR, middle row: SO, and bottow row: TO transitions. 


which is equal to 0 in the pure austenite and 1 in 
the pure/twinned martensite, 




= (3.1) 

r 

and specify conversion time t = tm when nm{tm) = 0.5. 

From Figure 3, we can see isothermal slow conversions 
and athermal fast conversions with incubation-delay tails 
for different material parameters Ai^Tc/Tq in TCR, SO, 
and TO transitions. Figure 3 also shows crossover from 
athermal to isothermal by fixing Ai and changing Tc/Tq^ 
and vice versa. Hence, we find the martensite classifica¬ 
tion is a matter of material parameters: the same model 
Hamiltonian can show both athermal or isothermal be- 
haviour, dependent on parameters. This is just as in the 
SR casM 

The athermal/isothermal/austenite regime diagrams 
are obtained in material parameters Tc/Tq, Ai and shown 
in Figure 4, that clearly depicts athermal martensites are 
more common than isothermaP. The simulations data 
matches well with the estimateJ^^ of theoretical bound¬ 
aries. Here, the criterion for athermal is tm = 10 MCS; 
isothermal/intermediate is tm = 1000 MCS; and austen¬ 


FIG. 4. Athermal/isothermal martensite and austenite phase 
regimes: Data of athermal, intermediate, and austenite be¬ 
haviour in a plot of 1 — ^ versus Ai, with Eq = 3 for TCR 
(top), SO (middle) and TO (bottom) transitions. Estimates of 
the theoretical boundaries are shown as solid lines. See text. 


ite, if there are no conversions even for holding time 
t = th. Again this is just as in the SR cas^l^. 

We will focus on the athermal regime. Figure 5 shows 
single-seed runs of n^(t) vs t after quenches to various 
Ar = T — < 0, below the transition temperature 

r 4 = t(T 4 ). At low temperatures, nm{t) rises rapidly 
to unity, but as transition is approached, shows incuba¬ 
tion behaviour. In the case of TCR transition, nm{t) 
rises to a smaller value, that incubates for longer times 
before it rises sharply to unity. In SO transition, n^(t) 
shows incubation followed by jerky steps before it rises 
to unity. In TO transition, n^(t) has longer incuba¬ 
tions before it sharply rises to unity. The transition 
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FIG. 5. Conversion incubation times: Martensite conversion fraction nm{t) versus time t, for stiffness Ai = 4 and various 
At = t — T4 < 0, showing the 50% conversion dehnition of tm for TCR (left), SO (middle) and TO (right) transitions. As 
athermal regime parameters have been chosen, flat incubation is seen for early times. 


is ’fuzzy’ and is operationally defined as where all 100 
runs give austenite. Hence, we defin^i^ mean conver¬ 
sion time = 1/ < >, where mean conversion rate 

< >=< 1/tm > is obtained by an arithmetic average 

over Nruns = 100 seeds. 

We henceforth focus on the athermal martensite pa¬ 
rameter regime to study the conversion-delays kinetics. 


IV. TEXTURAL ENERGIES PARAMETRIZED 
BY SURROGATE DROPLETS 

The transition is known to depend both on temper¬ 
ature and the size of martensitic seeds, as in the Pati- 
Cohen modeP. In early work, Pati and CoheiP have 
measured and modeled the conversion times in Ni-Mn al¬ 
loys and found that the isothermal slow conversions can 
change to athermal fast conversions, for fixed marten¬ 
site fraction, but with larger (and hence fewer) initial 
martensite seeds. This can be understood through the 
parametrization of textural droplet energies as in SR 
transitiorP^. At t = 0, the seeds of Ny variants are 
randomly sprinkled throughout the lattice. We find that 
the interaction tend to cancel leaving only self-interaction 
part Ai\U]/2 at each seed. So we have, 

^T[5(0)] 


the initial energy has a droplet-like form j3H{R{0)) = 
Co2^^i?c[l — (1 — i?(0)/i?c)^]. Here we define a length 
Rc{t) that is positive below a divergence temperature 
r = tl{Ai), 


Rc{t) 


- 2^2 

9L{r) + A,[U]/2- 


(4.2) 


As in the SR case, we define a scaled temperature variable 
?7(t) from the parametrization 


v{t) = -R{0)/Rc{t) = 


gLjr) + Ai[U]/2 

2 e 


(4.3) 


that will be used later, for i?(0) = 1. 

At t = 0, the initial seeds have a geometric meaning, 
and hence the pseudospin Hamiltonian H(S(0)) matches 
the droplet Hamiltonian H(R(0)), but for general t, these 
two terms no longer match term-by term. However, as 
in SR case, we define R(t) through H[S(t)]/H[S(0)] = 
H[R(t)]/H[R(0)]. The energy (ratio) for interacting vec¬ 
tor pseudospins is parametrized, by the energy (ratio) 
of a surrogate system of independent droplets. The ini¬ 
tially geometric R(0) evolves to an interacting-texture 
energetic parameter R(t), that can even go negative as 
the pseudospin energy goes negative. Thus 


p{t) = /3H{S{t))//3H{S{0)) 


+4Ms(f,of]. (4.1) 

Here, [U] ^ 0.5, is the Brillouin-zone average of U^^'(k) 
of (2.8) in TCR, SO, and TO transitions. For an ini¬ 
tial martensite fraction n^(0) = 0.02, we have Ny 
variants square seeds of sides R(0). The initial pseu¬ 
dospin seed energy is parametrized as /3H(S(0)) = 
Co[aLgLR{0)^ + aGe^R(0) + ac(A^[U]/2)R(0)^] with 
Co = (nm{0)NDo/2). For different sides R(0) = 1,2,3, 
we fit the coefficients aL^G,c term-by-term, finding again 
oll = OLG = Oic = 1, independent of seed size. Then 


= [1_(M_i)2]/[i_(«_i)2]_ (44) 

lijQ ILq 

The R(t) evolution is then once again 

R{t)/Rc = 1 + ay/l - pit)/Pc, (4.5) 

where pc = [1 ~ {{R{^)/Rc} — 1)^]”^? we take a = 
sign({R(0) / Rc} - 1). 

Figure 6 shows the evolutions of effective droplet en¬ 
ergy in a plot of R(t)/Rc(r) versus time. There are both 
rapid rises to final positive values and fiat-incubations as 
already seen in the martensite conversion fraction 
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FIG. 6. Trajectories : Scaled energy parameter R{t)/Rc ver¬ 
sus time t, showing flows are determined by initial R{0) / Rc(t) 
values. Note flat incubations of lower curves, of initial 1 > 
R(jd)/Rc(T) > 0.5, corresponding to T 2 < r < T 4 for TCR 
(left), SO (middle) and TO (right) transitions. 


which goes negative at later times. The flat-incubations 
are due to the inefficient searches for the rare channels 
to lower energies. The initial R(^)/Rc values determine 
the R(t) flows. 

As a consistency test of parametrization, Figure 7 
shows p(t)lpc versus R(t)/Rc indeed matches a parabola, 
for all t, and all Ai, and many starting values R(0)/Rc(r) 
in TCR, SO, and TO transitions. Flow directions of 
R(t) are indicated by arrows starting at R(0)/Rc for Re¬ 
gions 1 , 2 ,3,4, with asymptotic R(t) giving negative final 
martensitic energies, or zero (going to austenite). 



FIG. 7. Parametrization and crossover temperatures : Scaled 
pseudospin energy p(t)lpc versus R(t)/Rc(T), showing flows 
fall on a parabola as a test of parametrization. For R(0) = 
1 seeds, characteristic initial values R(0)/Rc are Rci~^ — 
2 ,Rc2~^ — 1,RcA~^ — 0.5 as marked. These correspond to 
temperatures ti,T 2 ,T 4 . For initial R(R) / Rc < 0.5, flows are 
to R = 0 austenite. 


i) Region 1: For initial R(0)/Rc(r) > 2, there are ex¬ 
plosive conversions to martensite, this determines a tem¬ 
perature r = Ti or 1/Rc(^i) = 2 or in scaled varibale 
77(r) = —2 with initial unit seeds R(0) = 1. 

ii) Region 2: For initial droplets 2 > R({))/Rc(t) > 1, 
the flows are again fast, this determines a temperature 
T = T 2 where 1 /Rc(^ 2 ) = 1 or r](r) = —1. 

hi) Region 3: For 0.5 > R(0)/Rc(r) or r](r) > —0.5, 


the initial droplets are flowing only to R = 0 austenite. 
But, for larger Ai, the droplets can still grow through 
searches up to R(0 )/Rc(t) ^ 0 or ? 7 (r) 0, that is well 

below the Landau transition temperature Tq. 

iv) Region f: For 0.5 < R(0) / Rc(r) < 1, the initial 
droplets immediately convert to a single variant droplet, 
that incubates for long times around R(t) 0 with zero 
energy H 0 (degenerate with austenite). This entrop- 
ically critical droplet searches for conversion pathways, 
and grows through jerky steps and autocatalytic twin¬ 
ning. The incubations occur for unit seeds up to a tem¬ 
perature r = T 4 or 1 /Rc(t 4 ) 0.5 or in scaled variable 

r](r) = —0.5. 

These critical values of the scaled variable r](T, Ai) are 
used in the scaled plot of Fig. 8 . 



FIG. 8. Textural crossover temperatures : Phase diagram 
Ro/Rc vs Ai. Theoretical boundaries (shown as dotted lines) 
in ti,T2,T4 are defined by Rc(ti,2,4) = Rci,c 2 ,c 4 of Fig 7 and 
symbols are data from simulations. See text. 


V. ATHERMAL REGIME CONVERSIONS 

In this Section, we find the conversion times and their 
distributions with a data collapse in terms of a scaled 
temperature-stiffness variable; and textural evolutions. 


A. Conversion time and their distributions 

The conversion times for different material parame¬ 
ters Ai,Eo and for TCR, SO, and TO transitions fall 
on a single hyperbola in logt^ versus At = r — T 4 , 
for a range of temperatures T 4 < r < ri as shown 
in main Fig 9. The same data is plotted in inset as 
l/ln(t^) versus |Ar|/|ri — T 4 I, that shows linearity on 
|Ar| goes to zero. The hyperbola and the linearity 
are showing the Vogel-Fulcher behaviomll^. Specifically, 
im = foexp[ 6 o|ri-T 4 |/|{r-r 4 }|], with to = 1.6, = IJ, 

for these data. The insensitivity of conversion times tm 
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FIG. 9. Singular divergence of conversion times: Plot of \ntm vs |At| = |t — T4I, for (a) TCR , (b) SO and (c) TO transitions 
with Tc/To = 0.9, Eq = 3,4, 5, 6 and Ai = 1,4,10. Solid line is In trn = In to + ^o|ti — T4|/| At|, with to = 1.6 and bo = 1.7. 


to energy scales Eq implies that the Vogel-Fulcher be¬ 
haviour at r 4 comes from divergence of entropy (rather 
than energy) barriers, tm ^ ^\^Sentr\ finding the rare 
channels^. The entropy barriers then vanish at ri, with 
a drop in conversion times. 

The main Figure 9 shows the temperature dependence 
of conversion time^^ for TCR, SO, and TO transitions. 
As in the scalar-OP SR transition, there are explosive 
conversions below a temperature r ^ ri (that is different 
for different transitions); and there are conversion delays 
for r > Ti, that rise at r r 2 , to diverge at a temperature 
r T 4 . The success conversion fraction versus Ar 
for various Ai = 1,4,10 and Eq = 3,4, 5, 6 with fixed 
Tc = 0.9 for TCR, SO, and TO transitions is shown in 
Supplemental Material. The fraction (frn is unity for r < 
r 2 and decreases for r > r 2 , to become zero at r T 4 . 
The insensitivity to different energy scales Eq is again a 
signature of entropy barrier^. 



0 0.1 0.2 0 0.1 0.2 0 0.1 0.2 


FIG. 10. Log-normal distribution of conversion-rates: Plot of 
P{rm) versus rm for different At. Data are from Fig.9, with 
Ai = 4,F^o = 3 for TGR (left), SO {middle) and TO {right) 
transitions. Solid lines are the log-normal distributions. 

Similar to the SR transition, we have calculated the 
arithmetic mean rate < >=< 1 /t^ > that deter¬ 
mines = 1 / < >, with 1 /t/i < < 1 in TCR, 

SO, and TO transitions. The variance in the rates is 
(Jr‘^ =< {vm— < Tm >)^ >• The probability densities 
P{Tm) versus for various Ar are shown in the Figure 
10, as histograms for different temperatures. For each 
his toy am of N^ist data points, the Scott optimized bin 
size^is used, of dvm = The histograms 


again narrow sharply, below r 2 , as in the delta-function¬ 
like peak on the right. Solid lines are Log-normal curves 
for calculated < Vm > and from the data. The Log- 

' m _ _I ^ 

normal distribution is a signature of rare events^. 



FIG. 11. Temperature-Time-Transformation phase diagram : 
Gonversion times tm versus scaled variable ? 7 (t) for TGR, SO, 
and TO transitions showing fast and delayed conversions for 
given parameters Ai = 4, Eq = 3, TcfTo = 0.9. 

Figure 11 shows TTT phase diagram in conversion 
times tm versus scaled variable 7^(r) = —1/Rc{r) for fixed 
Ai = 4^ Eq = 3 for TCR, SO, and TO transitions. The 
characteristic temperatures ti^ 2,4 are defined in terms of 
scaled variable as r = ri or 77(r) = —2 where tm — ^ 
MCS; r = r 2 or r]{T) = — 1 where tm — 10 MCS; and 
r = r 4 or ? 7 (r) = — 0 . 5 , where conversion times diverges. 


VI. EVOLUTION OF TEXTURES 

In the athermal parameter regime, after quenching into 
r 4 > r > r 2 , we monitor evolution of OP strain S tex- 















































FIG. 12. Evolution for the Triangle-centered rectangle (TCR) transition: First row: Snapshots of OP strain S = {S 2 ,S 3 ) for 
different times t on quenching to r = —2.7. See movie of this evolution. The color bar is in terms of variant label V. See text. 
Second row: Evolving stress distributions of P 2 {r),p 3 {r). Parameters: Ai = 4, Eq = 3. 



FIG. 13. Evolution for the Square-oblique (SO) transition: First row: Snapshots versus time t, of OP strain S = (*S' 2 ,*S' 3 ) on 
quenching to r = —1.1. The colour bar is in terms of variant label V. See movie. Second row: Evolving stress distributions 
P 2 {r),p 3 {f). Parameters: Ai = 4, Eo = 3. 


tures, local internal stresses (See Appendix.) and stress 
distributions to understand the conversion-incubations at 
microstructure level. The color bar in Figs 12, 13, 14, 15 
represents the OP strain S in terms of variant label V. 
In all the three transitions, ^ = 0 represents austenite 
S = 0 and V = l,2...Ny corresponds to Ny martensite 
variants with pseudospin vector values given in (2.14), 
and pictured in Fig 1. 

As shown in Figs. 12, 13, 14 (first row), after quench¬ 
ing the austenite with 2% martensite seeds into > 
r > r2, the seeds quickly form domain-wall Vapour’ of 
droplets of single variant (s), reminding Ostwald ripening. 
The droplets searches for the rare pathway channels to 
expand in the easy directions of the anisotropic potential. 
The expanded droplet then generates the other var iant 
by autocatalytic twinning as in Bales and Goodin^ZHH 
to form ’liquid’ of domain-walls, which orient to form 
domain-wall ’crystal’ at a later time. The jerkiness dur¬ 
ing conversion incubation is reflected i n wave number km 
(not shown) as steps with finite vadneJ^Eni ^^d also in 
(excess) thermodynamic quant it ie^^^^^, internal energy 


AU and entropy —TAS (not shown). 

In the second row of Figs. 12, 13, 14, the local stress 
distributions P 2 {'^)^P 3 {'k) are shown. At t = 0, the stress 
distributions are sharply peaked around zero with large 
values, which generate wings on both sides of the peak 
during autocatalytic twinning. In the final oriented state, 
only the wings remain, that correspond to the trapped 
stress values along the domain-walls (except in TO case, 
where P 3 {f) is sharply peaked around zero). 

The final ’equilibrium’ microstructures in the TTT 
phase diagram for TCR, SO, and TO transitions are 
shown in Figure 15 and are in good agreement with con¬ 
tinuous variable simulationand experimental^. 

With random initial seeds, there is a vibrating marten¬ 
site phase, that has bulk austenite in TCR, SO, and TO 
transitions as in Fig. 15 (a) , that could be equivalent of 
the chequerboard SR case tweed patterrlSl (and becomes 
less probable closer to r T4.). 

With 2% of martensite seeds, and for r > T4 or 7^(r) > 
—0.5, there is only uniform austenite. For r2 < r < r4 
or —1 < T] < —0.5, there are again austenite droplets 
but now appear as lines, in domain wall crystal (DWC) 
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FIG. 14. Evolution for the Triangle-oblique (TO) transition: First row: Snapshots of the OP strain S = {S 2 ,S 3 ) evolution in 
terms of V on quenching to r = — 1 . 1 . See movie. Seeond row: Stress distributions p 2 {T^,P 3 {r). Parameters: Ai = 4, = 3. 





(a) 


Bfia 


(b) (c) (d) (e) 



FIG. 15. Final ’equilibrium’ mierostruetures: The long time 
textures are shown, for transitions TGR, SO, and TO are in 
the first, second and third rows, respectively, (a) precursor 
vibrating phase (bulk-austenite) for r ^ T 4 ; (b) and (c) fan¬ 
like and Z-shaped twins (austenite appears as point or line 
densities) for T 2 ^ t > T 4 ; (d) star-like or martensite twins 
(no bound austenite) for n ^ t > T 2 ; (e) glass for r << n. 
The austenite is represented as zero in color bar. Parameters 
are typically Ai = 4, Tc = 0.9, = 1 and Eq = 3. 


in SO, TO cases, and Z-like stated in TOR case as in 
Fig. 15 (b). For ti < r < r2 or —2 < r]{r) < — 1, austenite 
droplets can appear as points at corners, in DWG (that 
have topological charges) in SO, TO cases; and also fan¬ 
like oriented stated in TOR case as shown in Fig. 15 (c). 
At low temperatures for r << ri or r]{r) « —2, the 
DWG or oriented twins can have vortex-like (or topolog¬ 
ical defects or charges) behaviour at multi-variant junc¬ 
tions in SCpl, TO cases, and partially oriented star-like 
state^I^ in TGR case, can compete with a frozen domain 
wall liquid or ’glass’ as shown in Fig. 15 (d), (e), and has 
no bound austenite. Hence, ri is austenite (local) spin- 
odal temperature. The microstructure (d) as shown in 
Fig. 15 (first row) for TGR transition is not fully relaxed 
even after th = 10^ MGS and could possibly take longer 


and longer times, to orient fully to a nested star as seen 
in continuous variable simulations and experimentJ^^^^I^. 


VII. CONVERSIONS WITHOUT THE 
COMPATIBILITY INTERACTION 


We now turn-off the compatibility term {Ai = 0) in the 
Hamiltonians for TGR, SO and TO transitions to under¬ 
stand the role of the power-law anisotropic potentials in 
the domain-wall conversion kinetics. 

Figure 16 showing the martensite fraction n^(t) ver¬ 
sus time t for different temperatures |Ar| = |r — r4|, 
where r4 is austenite transition temperature; and cover¬ 
sion times tm versus temperature At for different energy 
scales Eq for TGR, SO, and TO transitions with elastic 
stiffness constant Ai = 0. Here, rim has no flat regions 
or incubation as was seen for Ai = 4 in Fig.5. The final 
microstructure in Ai = 0 case is a slab-like martensite 
unlike oriented microstructures in Ai = 4 as in Fig. 15. 

Conversion times show a rise at ri, that almost re¬ 
main constant till r2 and then increase linearly to T4, 
above which there are no conversions found. There are 
no Vogel-Fulcher rises in conversion times as in Ai = 4 
(Fig.9), but there is a small Eq dependence, which could 
be now from energy barriers rather than entropy barriers. 

With the parametrization scale of (4.2) now given by 
Rc{Ai = 0) = 2 ^^/|^l|, the estimated transition temper¬ 
atures are ri = —5.5,—4.0,—4.0; r2 = —2.0,—1.4,—1.4; 
r4 = —0.4, —0.2, —0.2 are in good agreement with 

the simulation values of ri = —5.5,—4.0,—4.0; r2 = 
-2.0, -1.4, -1.4; r4 -0.8, -0.1, -0.1 for TGR, SO, and 
TO transitions respectively. 

Therefore, microstructures and conversion times in 
TGR, SO, and TO transitions with Ai = 0 are clearly 
different from the non-zero Ai case. Thus, the power- 
law anisotropic potentials in ferroelastic transitions are 
important in understanding orientations and kinetics of 
domain-walls. 
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FIG. 16. Micro structures and conversion times with Ai = 0; 
The left column shows martensite fraction nm{t) versus t for 
|At| = 0.2,0.3,0.5 and Eq = 3, with the conversion time 
tm is marked for TCR (top row), SO {middle row), and TO 
{bottom row). The right column shows the conversion times 
im versus At for Eq — 3,4,5. Compare to Fig 5 and Fig. 9, 
respectively. 


VIII. SUMMARY AND FURTHER WORK 


Systematic temperature-quench MC simulations with¬ 
out extrinsic disorder are carried on the strain pseudospin 
clock-zero model Hamiltonians, with vector-order param¬ 
eter {Nqp = 2) and Ny -\-l{= 4, 5, 7) strain-pseudospins, 
that correspond to triangle-to-centred rectangle, square- 
to-oblique, and triangle-to-oblique transitions to get in¬ 
sights into conversion-incubation kinetics. The results 
are similar to the SR cas^^, that are just seen to be 
generic. The simulation results are as follows: 

(1) The microstructures of discrete-strain pseudospins 
in the Temperature-Time-Transformat ion phase diagram 
are in good agreement with continuous-variable simula¬ 
tions and experiment. 

(2) On quenching, for different material parameters 
Tc/To,^‘^, Ai, Eq, martensite fraction nm{t) can have 
slow isothermal and fast athermal conversions. The con¬ 
version times tm can transform from rapid athermal to 
slow isothermal or vice versa on changing the material 
parameters; and athermal/isothermal/austenite regime 
diagrams are obtained in material parameters. 

(3) Focusing on the athermal parameter regime, we 


find rapid conversions below a spinodal like temperature 
and incubation delay above it, as in the experiment. The 
conversion-delay times have Vogel-Fulcher divergences, 
which are insensitive to Hamiltonian energy scales Eq 
and conversion rates have Log-normal distributions, as 
in scalar-OP SR transition, from entropy barriers. 

(4) The Temperature-Time-Transformation diagram in 
the athermal parameter regime has crossover tempera¬ 
tures and are understood through parametrization of do¬ 
main wall textures by surrogate droplet energies. 

(5) During conversion incubation (t^), evolutions of 
microstructures, stress distributions and domain-wall 
thermodynamics are monitored. The initial marten¬ 
site seeds in the austenite at t = 0 disappear to form 
a domain-wall vapour of single variant droplet (s), that 
incubates before generating Ny variants, one after the 
other, by autocatalytic twinning to convert to domain- 
wall liquid. The wandering domain-walls then orient 
later to a domain-wall crystal. During incubation, stress 
distributions remain sharply peaked, and there are finite 
steps in (excess) internal energy, (excess) entropy. 

(6) On switching off the power-law anisotropic poten¬ 
tials, we find no incubations in conversions, no Vogel- 
Fulcher divergences and the microstructure is multi-slab 
martensite. 

Systematic experiments on athermal martensites can 
look for martensite formation and growth during conver¬ 
sion incubation and their divergences as well as distribu¬ 
tions close to the transition. 

Monte Carlo simulations presented in this paper are 
on 2D strain-pseudospin models for ferroelastic transi¬ 
tions with vector-OP. We also find similar conversion 
incubation-delays in 3D strain-pseudospin models for 
tetragonal-to-orthorhombic {Ny = 2, Nqp = 1), cubic- 
to-tetragonal {Ny = 3, Nqp = 2), cubic-to-orthorhombic 
{Ny = 6 , Nqp = 2), and cubic-to-trigonal {Ny = 
4, Nop = 3) ferroelastic transitions^. 
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APPENDIX: INTERNAL STRESSES 

The local internal stresses, P2(^) = dF/de 2 {r)^ and 
P3(r) = dF/des{r) for TCR transition are obtained as, 

P 2 {r) = e[4e^S'| — 6e5'| + 25'2(2e^5'| + r) -k 6e5'| 


+2^^A^5'2 + Ai{U 22 S 2 + UqqSq)], (A.la) 


P3(^) = e[4e^5'| + 2S'3(2e^5'| + 6e5'2 + r) + 2^^A^S'3 


+Ai(U23^2 + U33^3)]. {A.lb) 




























11 


The local internal stresses for SO transition are, 
P2{r) = e[2S'2(3S'2 - 4S'| + r) + 2^^A^S'2 

+Ai{U22S2+Us 3S3)], (A.2a) 

Pz{r) = e[2S'2(3S'| - 4S'| + r) + 2^^A^S'3 

+Ai([/2352 +C/ 3353 )] • {A.2b) 

The local internal stresses for TO transition are, 
P2{r) = e[252(35| - 45| + r) + 2C^A^52 


+Ai(C/2252+C/ 3353)], (A.3a) 


Psi''’) = £[252(353 - 45| + t) + 2^^A^53 


+Ai(C/23 52 +C/ 3353)]. (A36) 
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